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^! . ABSTRACT 

A Suzaku observation of the nucleus of the radio-loud AGN Centaurus A in 2005 has yielded a broad- 
. band spectrum spanning 0.3 to 250 keV. The net exposure times after screening were: 70 ks per X-ray 

Imaging Spectrometer (XIS) camera, 60.8 ks for the Hard X-ray Detector (HXD) PIN, and 17.1 ks for 
C**"*- , the HXD GSO. The hard X-rays are fit by two power laws of the same slope, absorbed by columns 

of 1.5 and 7 x 10 23 cm -2 respectively. The spectrum is consistent with previous suggestions that the 
power-law components are X-ray emission from the sub-pc VLBI jet and from Bondi accretion at the 
core, but it is also consistent with a partial-covering interpretation. The soft band is dominated by 
thermal emission from the diffuse plasma and is fit well by a two-temperature VAPEC model, plus a third 
power-law component to account for scattered nuclear emission, jet emission, and emission from X-ray 
Binaries and other point sources. Narrow fluorescent emission lines from Fe, Si, S, Ar, Ca and Ni are 
detected. The Fe Ka line width yields a 200 light-day lower limit on the distance from the black hole 
to the line-emitting gas. Fe, Ca, and S K-shell absorption edges are detected. Elemental abundances 
are constrained via absorption edge depths and strengths of the fluorescent and diffuse plasma emission 
lines. The high metallicity ([Fe/H]=+0.1) of the circumnuclear material suggests that it could not have 
originated in the relatively metal-poor outer halo unless enrichment by local star formation has occurred. 
Relative abundances are consistent with enrichment from Type II and la supernovae. 

Subject headings: galaxies: active — X-rays: galaxies — galaxies: individual (NGC 5128) 



1. INTRODUCTION 

The radio-loud active galactic nucleus (AGN) Centau- 
rus A (NGC 5128) is one of the most extensively studied 
AGNs at all wave bands, thanks to its proximity (distance 
of 3.8 ± 0.4 Mpc, Rejkuba 2004; V = 1 kpc) and bright- 
ness (2-10 keV flux typically ~2 xl0~ 10 erg cm~ 2 s^ 1 ). 
Cen A features the nearest AGN jet and is considered a 
prototypical FR I radio galaxy. Its radio structure con- 
sists of extended kpc-scale outer, middle, and inner radio 



lobes, fed by a one-sided kpc-scale jet. The jet is at an 
angle of ~60° to the line of sight, leading to its descrip- 
tion as a 'misaligned' BL Lac object (e.g., Bailey et al. 
1986). VLBI observations have also revealed a sub-pc jet 
and unresolved (~ O.lpc) core (Tingay et al. 1998; see also 
Israel 1998 for a review). 

The host galaxy of Cen A is a giant elliptical with a pro- 
nounced optical dust lane that obscures the inner few kpc 
and is believed to be an edge-on disk structure (Quillen et 
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al. 1992). The inner dust disk, with associated young stars, 
H II regions, and other ionized gas, is the site of vigorous 
star formation, e.g., observed along the edges of the in- 
ner disk (e.g., Ebneter & Balick 1983, Dufour et al. 1979). 
The dust disk is interpreted as the result of a merger be- 
tween a large elliptical and a gas-rich disk galaxy at least 
a Gyr ago (e.g., Israel 1998). Infrared studies are neces- 
sary to reveal the compact, pc-scale core (e.g., Karovska 
et al. 2003). NIR studies (e.g, Schreier et al. 1998) have 
revealed a circumnuclear disk tens-hundreds of pc wide, 
surrounding the central black hole. Stellar kinematic stud- 
ies indicate a black hole mass near 2xl0 8 M (Silge et al. 
2005) 

Previous studies of the nuclear X-ray emission from —4 
keV to several hundred keV (e.g., Rothschild et al. 1999) 
have established the presence of a non-thermal, power- 
law continuum whose origin is not certain. It could be a 
signature of accretion, or it could be associated with jet 
emission processes, e.g., synchrotron or inverse Compton 
emission from the sub-pc VLBI jet. 

Below —4 keV, the X-ray continuum emission under- 
goes moderately heavily absorption, although the nature 
of the X-ray obscuring material is not certain; i.e., it could 
have a sky-covering fraction of 1 as seen from the central 
X-ray source or it could be in the form of a dusty torus. 
Rothschild et al. (2006) have noted —50% variations in 
the absorbing column over —20 years, and suggested that 
the nucleus is seen through the edges of a warped, ro- 
tating disk. The strong X-ray absorption makes Cen A 
relatively unique: the vast majority of FR I radio galax- 
ies possess relatively weak X-ray absorbing columns (< 
i()20-2i cm -2 ). Heavy X-ray absorption, meanwhile, is 
usually seen only in FR II radio galaxies, although it is 
not yet clear whether the dichotomy in X-ray absorption 
properties has the same origin as the FR I /ll dichotomy, 
e.g., differing accretion modes (Donato et al. 2004, Bal- 
maverde et al. 2006, Evans et al. 2006). 

Previous X-ray spectral fits to the nucleus of Cen A 
center on one or more power laws, with varying degrees 
of absorption. For instance, using ROSAT and ASCA, 
Turner et al. (1997) suggested that one plausible descrip- 
tion of the nuclear emission is a partial-covering model, 
wherein the nucleus is seen through three layers of ab- 
sorption: 40%, 59%, and 1% of the nuclear emission is 
obscured by columns of 4, 1 and 0.01 x 10 23 cm~ 2 , re- 
spectively. Using Chandra-EETGS and XMM-Newton ob- 
servations, Evans et al. (2004) suggested a different in- 
terpretation. They fit the hard X-ray continuum spec- 
trum with a heavily-absorbed power law, thought to be 
associated with Bondi accretion at the core, plus a rela- 
tively less-absorbed power law thought to be associated 
with the sub-pc VLBI jet, although they did not rule out 
the partial-covering interpretation. 

There is a prominent Fe Ka emission line at 6.4 keV, 
first noted by Mushotzky et al. (1978). It is narrow (e.g., 
FWHM = 2200±900 km s" 1 , Evans et al. 2004), and its 
flux is historically nearly constant despite long-term con- 
tinuum variations spanning a factor of —5 over 20 years 
(Rothschild et al. 1999), implying an origin for the line 
that is distant from the origin of the variable continuum. 

The soft X-ray emission within —6 kpc of the nucleus 
of Cen A consists of both diffuse and point-like emission. 



There is diffuse emission associated with the kpc-scale jet 
extending from < 60 pc to the NE radio lobe (e.g., Kraft 
et al. 2002), as well as from thermal gas that surrounds 
the nucleus out to a radius of roughly 6 kpc. Point-like 
emission is associated with several dozen knots in the kpc- 
scale jet (e.g., Kataoka et al. 2006) as well as from a popu- 
lation of X-ray Binaries (e.g., Kraft et al. 2000). A soft X- 
ray spectrum covering the inner few kpc is thus expected 
to contain both power-law continuum emission from the 
jet and point-like sources and line- like emission associated 
with the thermal gas. For instance, Turner et al. (1997) 
modeled the soft X-rays with a 0.6 keV thermal plasma 
component, a 5 keV component associated with emission 
from a population of X-ray Binaries, and a power-law jet 
component. 

In this paper, we report on an observation of Cen A 
made with the Suzaku observatory in 2005 August. The 
combination of the X-ray Imaging Spectrometer (XIS) 
CCD and the Hard X-ray Detector (HXD) instruments 
have yielded a broadband spectrum covering 0.3 to 250 
keV, allowing us to deconvolve the various broadband 
emitting and absorbing components in this object. Fur- 
thermore, the exceptional response of the XIS CCD and 
high signal-to-noise ratio of this observation allow us to 
study narrow emission lines in great detail. This observa- 
tion has also yielded the highest quality soft X-ray spec- 
trum of Cen A obtained so far. §2 gives a brief overview of 
the Suzaku observatory, and describes the observation and 
data reduction. §3 describes the spectral fits. The results 
are discussed in §4, and a brief summary is given in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

The nucleus of Cen A was observed by Suzaku from 2005 
August 19 at 03:30 UT until August 20 at 09:50 UT, and 
in fact was the first light target for the HXD. Suzaku was 
launched 2005 July 10 into a low-Earth orbit. It has four 
X-ray telescopes (XRTs; Serlemitsos et al. 2007), each with 
a spatial resolution of 2' (half-power diameter) . The XRTs 
focus X-rays onto four X-ray Imaging Spectrometer (XIS; 
Koyama et al. 2007) CCDs, which are sensitive to 0.2- 
12 keV X-rays on a 18' by 18' field of view, contain 1024 
by 1024 pixel rows each, and feature an energy resolu- 
tion of -140 eV at 6 keV. Three CCDs (XIS0, 2, and 
3) are front-illuminated (FI), the fourth (XIS1) is back- 
illuminated (BI) and features an enhanced soft X-ray re- 
sponse. The XRT/XIS combination yields effective areas 
per detector of roughly 330 cm 2 (FI) or 370 cm 2 (BI) at 
1.5 keV, and 160 cm 2 (FI) or 110 cm 2 (BI) at 8 keV. Each 
XIS is equipped with two 55 Fe calibration sources that 
produce fluorescent Mn Ka and Kf3 lines and are located 
at the CCD corners. Suzaku also features a non-imaging, 
collimated Hard X-ray Detector (HXD; Takahashi et al. 
2007); its two detectors, PIN and GSO, combine to yield 
sensitivity from ~10 to —700 keV. Further details of the 
Suzaku observatory are given in Mitsuda et al. (2007). 

2.1. XIS reduction 

The XIS data used in this paper were version 0.7 of 
the screened data (Fujimoto et al. 2007) provided by the 
Suzaku team. The screening is based on the following crite- 
ria: grade 0, 2, 3, 4, and 6 events were used, the CLEANSIS 
script was used to remove hot or flickering pixels, data col- 
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lected within 256 s of passage through the South Atlantic 
Anomaly (SAA) were discarded, and data were selected 
to be 5° in elevation above the Earth rim (20° above the 
day-Earth rim). The XIS FI CCDs were in 2x2, 3x3 and 
5x5 editmodes, for a net exposure time after screening of 
71.5 (XISO), 70.3 (XIS2) and 68.6 (XIS3) ks. XIS1 was in 
3x3 and 5x5 modes, for a net exposure of 68.6 ks. 

The source was observed at the nominal center position 
of the XIS. The XIS was kept in Normal Mode, but to 
reduce the risk of photon-pileup, the 1/4 Window Option 
was used: each CCD recorded 256 by 1024 pixel rows, 
with a readout time of 2 s, as opposed to 8 s when the full 
1024 by 1024 window is used. For each XIS, we extracted 
a 2' radius centered on the source. Spectra were binned 
to a minimum of 50 counts bin -1 to allow use of the x 2 
statistic. 

The count rates observed over the full XIS band were 5.0 
ct s- 1 (average of XISO, 2 and 3) or 4.3 ct s^ 1 (XIS BI). 
These count rates are far below the limit for photon-pileup 
in 1/4 Window Mode, ~ 12.5 ct s _1 . Furthermore, ex- 
tracting over annular regions with inner and outer radii of 

0. 4' and 2', respectively, yielded spectra identical in shape 
to those extracted over 2' circular regions, indicating that 
pile-up was negligible. 

Because the soft diffuse emission in Cen A covers nearly 
the entire read-out area of the CCD, we used observations 
of the north ecliptic pole (NEP) to extract a background 
spectrum. Suzaku observed the NEP on 2005 September 
2 to 4. There was a brief, sharp flare during the NEP ob- 
servation, consistent with evidence for solar wind charge 
exchange emission (see Fujimoto et al. 2007). Removing 
that flare and applying the same screening criteria as above 
yielded a net exposure time of 93.9 (XISO), 108.5 (XIS1), 
93.3 (XIS2), and 96.6 (XIS3) ks. We extracted data over a 
2' circle. The NEP is known to contain soft X-ray emission 
lines (see e.g., Fujimoto et al. 2007), but they are much 
fainter than those of Cen A (see §3). The O VII and O 
VIII lines in Cen A are a factor of roughly 10 (20) higher 
compared to the NEP as seen in the XIS BI (XIS FIs). 
Lines near 0.8-0.9 keV such as Fe L XVII and Ne IX in 
Cen A are typically —60 (100) times brighter in Cen A 
compared to the NEP as seen in the XIS BI (XIS FIs). 

Average 0.5-2.0 keV count rates for the NEP observa- 
tion were 0.004, 0.009, 0.004, and 0.003 count s" 1 for XISO, 

1, 2, and 3, respectively; average 2-10 keV count rates were 
0.005, 0.012, 0.005, and 0.004 count s" 1 , respectively. 

Response matrices and ancillary response files (ARFs) 
were generated for each XIS independently using xissim- 
RMFGEN and XISSIMARFGEN version 2006-10-26 (Ishisaki 
et al. 2007). The ARF generator takes into account the 
level of hydrocarbon contamination on the optical block- 
ing filter. However, the Cen A observation occurred very 
early in the Suzaku mission, and the level of contamina- 
tion was quite low: we estimate a carbon column density 
of only -0.1, 0.2, 0.3 and 0.5 xlO 18 cm~ 2 for XISO, 1, 2 
and 3, respectively. 

A consequence of observing with the 1/4 Window Op- 
tion is that there are no 55 Fe calibration source lines on 
the CCD. This leads to a higher systematic uncertainty in 



the energy scale compared to normal mode (where it is a 
few tenths of a percent at most) . To estimate the instru- 
ment resolution, we used the calibration source lines from 
a 77 ks observation of MCG-6-30-15, observed in Normal 
Mode with no window option immediately prior to Cen A. 
We fit the calibration source spectra with three Gaussians. 
Two Gaussians were for the Mn Ka doublet (expected en- 
ergies 5.899 keV and 5.888 keV), with energy centroids 
fixed to be 11 eV apart, and the higher energy line flux set 
to twice that of the lower energy one. The third Gaussian 
was used to model the K/3 line, expected at 6.490 keV. We 
found the average of all the calibration line widths a to be 
9^2 eV; residual width may be due e.g., to imperfect CTI 
correction. 

The positional accuracy of Suzaku is ~1' at present, as 
the spacecraft has been known to exhibit small attitude 
variations ("wobble"). We checked to make sure that the 
effect of these attitude variations was not significant, given 
that the 1/4 Window Option was used throughout the ob- 
servation. We generated light curves by extracting over 
radii of 1.5' and 1.0' and looked for variations on the or- 
bital timescale, but found nothing significant. 

Average 0.5-2.0 keV net source count rates were 0.24, 
0.35, 0.25, and 0.25 count s" 1 for XISO, 1, 2, and 3, re- 
spectively. Average 2-10 keV net source count rates were 
4.7, 4.4, 4.5, and 4.5 count s" 1 for XISO, 1, 2, and 3, re- 
spectively. Figure 1 shows 2-10 keV light curves for each 
XIS. During the observation, the 2-10 keV source flux in- 
creased by only 10%, with fractional variability amplitudes 
-Fvar (see Vaughan et al. 2003 for definition) of 3.7 ± 0.2 
%, 3.4 ± 0.2 %, 3.1 ± 0.2 %, and 3.3 ± 0.2 % for XISO, 1, 
2, and 3, respectively. The 0.5-2.0 keV flux (determined 
from model fits; see §3) was 2.33xl0~ 12 erg cm~ 2 s _1 . 
The 2-10 keV flux was 2.12xl0~ 10 erg cm" 2 s _1 . 



2.2. HXD reduction 

The PIN and GSO source spectra were extracted from 
cleaned version 1.2 (pre-1.2-rl) HXD event files provided 
by the HXD instrument team. We first discuss the PIN ex- 
traction. Data were selected according to the following cri- 
teria: at least 500 s since SAA passage, COR > 8 GV, and 
day- and night-Earth elevation angles each >5°. Instru- 
mental (non- X-ray) background spectra for the PIN were 
provided by the HXD Team ( "background D" model) gen- 
erated from a time-dependent model(3 The model utilized 
the count rate of upper discriminators as the measure of 
cosmic ray flux that passed through the silicon PIN diode 
and yielded background spectra based on a database of non 
X-ray background observations with the PIN (Fukazawa 
et al. 2007). The systematic uncertainty of the PIN is ex- 
pected to be <5%. However, we note that Cen A is one 
of the brightest AGNs in the PIN energy range, and the 
effect of background subtraction error on the net source 
spectrum is relatively small. Both the source and back- 
ground spectra were generated with identical good-time 
intervals, and the exposures were corrected for instrument 
dead time (a ~5% effect). This yielded a good-time expo- 
sure of 60.8 ks. Data < 12 keV were discarded due to noise 
contamination near the lower threshold of the PIN diode. 
Data above 76 keV were also discarded: the gain above 

13 Other recently-published Suzaku results have made use of the "background A" model. However, observations occurring before 2005 Septem- 
ber 2 had a hit-pattern width set to a shorter value compared to observations occurring after this date. Because applying the "background A" 
model would yield an underestimate of the true non- X-ray PIN background, we instead use the "background D" model. 
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an internal Bi Ka calibration line at 76 keV is not well- 
defined, although there are GSO data covering these ener- 
gies. Further details of the HXD in-orbit performance are 
given in Kokubun et al. (2007) . To model the contribution 
to the total background from the cosmic x-ray background 
(CXB), the spectrum of the form 9. Ox 10" 9 (E/3keV)-°- 29 
exp(-E/40keV) erg cm" 2 s" 1 sr" 1 keV" 1 (Gruber et al. 
1999) was used; the contribution in the 12-76 keV band 
was l.lxlO -11 erg cm -2 s" 1 . 

The spectrum was binned to a minimum of 
400 count bin -1 . We used the response hie 
ae_hxd_pinxinom_20060814.rsp. The 12-76 keV net source 
flux and count rate were 7.3 xlO -10 erg cm -2 s _1 and 1.21 
count s" 1 . The total (X-ray plus particle) background 12- 
76 keV flux and count rate were 6.9xl0 -10 erg cm" 2 s" 1 
and 0.59 count s" 1 . The orbitally-binned light curve is 
shown in Figure 1. The light curve generally increases 
in flux by ~10-20% over the observation, matching the 
general trend shown by the XIS light curve, although 
shorter-timescale differences between the XIS and PIN 
light curves might be attributed to uncertainty associated 
with the PIN background. Figure 2 shows the net source, 
background, and total (source + background) spectra. 
The source spectrum is always at least 30% of the total 
up to -50 keV. 

The HXD GSO data were reduced in a manner similar 
to the HXD PIN data. A key contributor to the GSO par- 
ticle background is activation lines, e.g., delayed emission 
from radioactive isotopes induced inside the detector due 
to interaction with SAA particles. To minimize such emis- 
sion, data taken within 6000 s of SAA passage were dis- 
carded; the background is the most reliable and the least 
variable during non-SAA orbits. This yielded a good-time 
exposure of 17.1 ks. The 'background D' model hies pro- 
vided by the HXD Team were used (Fukazawa et al. 2007). 
The CXB is expected to contribute insignificantly to the 
total GSO background and was ignored. Source and back- 
ground spectra were both binned with GRPPHA following 
recommendations from the HXD TeaixQ 

The source was detected out to 250 keV; the GSO is 
sensitive down to roughly 45 keV. The 45-250 keV net 
source flux and count rate were 7.2xl0 -10 erg cm" 2 s" 1 
and 0.83 count s -1 . The 45-250 keV background flux and 
count rate were l.OxlO -8 erg cm" 2 s" 1 and 9.1 count 
s" 1 . Figure 1 shows the GSO light curve. Figure 2 
shows the net source, background, and total (source + 
background) spectra. The source spectrum is always at 
least 5% of the total below 200 keV. The response hie 
ae_hxd_gsoxinom_20060321.rmf was used. 

Below 100 keV, the GSO held of view is 34'x34' FWHM, 
the same as the PIN. Above 100 keV, the held of view in- 
creases with photon energy to a maximum of 4.5° x4. 5° 
FWHM and there is the possibility of source confusion. 
The nearest possible contaminating source is the blazar 
MS 1312.1-4221 = GRO J1312-42, located ~2° west of 
Cen A. According to Steinle et al. (1998), the flux of 
MS 1312.1-4221 in the Compton Gamma-Ray Observatory 
(CGRO) OSSE band (50-4000 keV) is likely less than a 

14 group 24 25 25 26 2 27 28 2 29 31 3 32 35 4 36 38 3; group 39 42 4 
92 100 9; group 101 110 10 111 121 11 122 134 13 135 147 13 148 162 15 



tenth of Cen A's flux. MS 1312.1-4221 is much fainter at 
other X-ray bands as well: the 1 keV intensity is usually 
only 1-10% of the historically observed Cen A intensities 
(Kinzer et al. 1995). We henceforth assume that the con- 
tamination from MS 1312.1-4221 is negligible in the GSO 
band. 

3. SPECTRAL ANALYSIS 

Given the 2' half-power diameter of the XRT, the spec- 
trum will contain blended contributions from multiple 
sources. Previous studies (e.g., Evans et al. 2004) have 
shown the hard X-ray^3 to be dominated by two power 
laws; both hard X-ray components are constrained to lie 
within a few tenths of an arcsecond (~5 pc) of the position 
of the unresolved radio core by Chandra ACIS (Kraft et 
al. 2000). Contributions to the soft X-rays include thermal 
emission from the diffuse plasma plus diffuse emission from 
the kpc-scale jet, with contributions from a few of the in- 
nermost knots in the kpc-scale jet and other point sources 
resolved by Chandra ACIS (Kraft et al. 2000; Kataoka et 
al. 2006). Each knot, however, is extremely faint in the 
soft X-rays (< 6 x 10" 13 erg cm" 2 s , Feigelson et al. 
1981). 

The four XIS spectra were fit separately; 0.4-11.5 keV 
FI data and 0.3-11.5 keV BI data were included. We 
found it reasonable to keep the photon indices tied for 
all four XlSes. while allowing the relative normalizations 
for XIS0 and 1 to each vary relative to that for XIS2 and 
3, which were tied together (XIS3/XIS2 was always con- 
sistent with, and therefore hxed at, 1.00). XIS0/XIS2 was 
typically 1.06; XIS1/XIS2 was typically 0.96 (in our best- 
fit model, XIS0/XIS2 = 1.065±0.004 and XIS1/XIS2 = 
0.951±0.003). We ignored 1.80-1.87 keV in the XIS FI 
spectra and 1.75-1.87 keV in the XIS BI spectrum due 
to uncertainties in calibration associated with the instru- 
mental Si K edge. This means we cannot directly study 
the intrinsic neutral Si K edge in the absorbing mate- 
rial or any highly-ionized Si K lines expected due to the 
thermal plasma emission. The PIN/XIS2 and GSO/XIS2 
normalizations were left free, but were usually close to 
1.09 and 0.9-1.0, respectively, although results were not 
strongly dependent on these factors (in our best-fit model, 
PIN/XIS2 = 1.055±0.022 and GSO/XIS2 = 0.98±0.06). 
Preliminary results on the relative GSO/PIN normaliza- 
tion using' the Crab suggest GSO/PIN = 0.84 (T. Ya- 
masaki et al. 2007, in preparation). Fixing GSO/PIN at 
0.84 in the Cen A hts usually resulted in % 2 increasing 
by <5; the photon index did not change significantly. We 
left the GSO/PIN normalization free to attain the lowest 
X 2 possible. All errors on one interesting parameter cor- 
respond to A\ 2 — 2.71 (with all relative normalizations 
except XIS3/XIS2 left free). Observed spectral features in 
Cen A are redshifted by z = 0.001825 (Graham 1978), the 
value corresponding to Cen A's recessional velocity. To 
calculate luminosities, a redshift of 0.0009, corresponding 
to a luminosity distance of 3.8 Mpc (assuming Hq = 70 km 
s" 1 Mpc" 1 ) was used. Galactic absorption of 8.6xl0 20 
cm" 2 was included. The abundances of Lodders (2003) 

43 46 4 47 51 5 52 56 5 57 62 6; group 63 68 6 69 75 7 76 83 8 84 91 8 
; group 163 178 16 179 196 18 197 216 20 217 238 22 239 262 24; group 
43 466 511 46 



263 288 26 289 317 29 318 349 32 350 384 35 385 422 38; group 423 465 
15 We henceforth denote "hard X-ray" emission as any intrinsic emission component above 2—3 keV, as opposed to describing emission detected 
solely with the HXD. 
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were used. Absolute abundances relative to solar values 
are denoted by Zp e ; a value of 1 denotes solar abundance. 

As shown below, the best model fit includes five over- 
lapping broadband components, and changes in how one 
is fit can potentially cause changes in other components. 
Our strategy in finding the best-fit broadband model was 
to fit the hard band first, then extend it to the soft band. 
Relatively more narrow features (lines and edges) were ex- 
plored later. 

3.1. Broadband fit 

We first fit the >3 keV emission with an absorbed power 
law (henceforth denoted PL1) using zvphabs(cutoffpl) 
(Model 1). The power-law cutoff was kept fixed at 400 
keV; abundances in the absorber were initially kept fixed 
at solar. Residuals are plotted in Figure 3. The slight 
hard excesses above ~9 keV in the X1S spectra are asso- 
ciated with calibration uncertainties known at the time of 
this writing. Residuals in the GSO spectrum near 150-200 
keV are likely associated with uncertainties in estimating 
the strengths of activation lines in the GSO background 
(lines at 153 and 196 keV are due to 153 Gd and 151m Eu, 
respectively); see Kokubun et al. (2007) for further details. 
GSO residuals near 65 keV are likely associated with cal- 
ibration uncertainties present at the time of this writing. 
Large residuals near the expected Fe Ka line were appar- 
ent. We then added Fe Ka and K/3 lines; the K(3 energy 
was fixed at 7.056 keV, its intensity was fixed at 0.13 times 
that of Ka, and the widths of the two lines were tied. The 
X 2 /dof dropped from 10700/7135 to 7607/7132, but the 
residuals in the XIS band suggested the fit could be im- 
proved if we added a second, fainter absorbed power law 
(PL2; Model 2). The photon index of the second power 
law, T 2 , was kept tied to that of PL1, Ti, for simplicity, 
but the absorbing column densities Ahi and A^H2 were 
not. x 2 dropped by over 140. PL1, with a 1 keV normal- 
ization ~4 times that of PL2, was absorbed by a column 
JVm near l.OxlO 23 cm~ 2 ; A H2 was near 5.0xl0 23 cm~ 2 . 
The photon index was 1.7T3±2!oii- Residuals to Model 2 
are plotted in Figure 3. 

We then included the soft band data, down to 0.3 keV. 
Several soft X-ray emission lines are prominent, and are 
especially well constrained with XIS1, as illustrated in Fig- 
ure 4. Simply adding a third power law (Model 3) to the 
soft X-rays thus did not yield an acceptable broadband fit. 
Table 1 shows results for the lines when each is fit with 
a Gaussian; we identify them as originating in O VII, O 
VIII, a blend of Fe XVII 3s-2p lines (the 3 Pi/ 3 P 2 dou- 
blet and the 1 P 1 line), a blend of Fe XVII 3d-2p lines 
{ l P 1 and 3 L>i), Ne IX, Ne X, and Mg XI. We can use the 
lines as a check of the absolute energy scale of the XIS. 
Assuming that in each He-like triplet, the resonant line is 
a factor of 2 in intensity greater than the forbidden line 
(roughly consistent with the plasma temperatures in our 
best-fit model, below), the measured energy centroids are 
within <16 cV of the expected energies. In addition, the 
Fe Ka centroid is 6±3 eV from its expected energy. 

The emission lines could be due to gas that is photo- 
or collisionally ionized. A key to distinguishing between 
the resulting spectra is that collisional ionization usually 
leads to the presence of very strong Fe L emission, such 
as a narrow Fe XVII emission line near 0.83 keV and a 



Fe L emission "bump" near 0.7-0.9 keV. We modeled the 
soft emission using XSTAR tables appropriate for emis- 
sion from photoionized gas, assuming either one or two 
zones of photoionized gas. Fitting this model always left 
strong residuals due to improperly modeled Fe L emission, 
even when the Fe abundance was left as a free parameter 
(in which case values commonly pegged at an upper limit 
of 10). The best-fit models assuming either one or two 
zones had x 2 /dof values of 11600/9159 and 11000/9155, 
respectively. 

Attempting to model the soft band emission with 
a single-temperature vapec component, with all abun- 
dances fixed at solar, also did not yield an acceptable fit. 
In Model 4, we tried two vapec components, henceforth 
denoted VAPEC1 and VAPEC2. The best-fit model had 
temperatures of k^T ~ 0.2 keV, to model the He- and H- 
like O lines, and k-gT ~ 0.6 keV, to model the Fe L and 
He- and H-like Ne lines. However, x 2 I dof was still high, 
10395/9160, and there were still large residuals <0.7 keV. 

It is plausible that some optically-thin circumnuclear 
material may scatter some of the hard X-ray continuum; a 
scattered power-law component was included in the best- 
fit model of Turner et al. (1997), for instance. In addition 
to the two vapec components, we added a third power law 
(PL3), with photon index T 3 tied to Ti, assuming that PL3 
undergoes no absorption in excess of the Galactic column. 
Here x 2 /dof was high: 10158/9161. Untying T 3 improved 
the fit considerably: x 2 /dof fell to 9806.8/9158; T 3 was 
1.16±8;^ (Model 5). 

However, it is plausible that PL3 does undergo some ab- 
sorption in the host galaxy, so we added a zvphabs com- 
ponent to PL3, and refit with r 3 tied to Ti. x 2 /dof fell 
to 9798.9/9160. The column density of the third absorber, 
A H3 , was 0.13±0.06 xlO 22 cm" 2 . 'The 1 keV normaliza- 
tion of PL3 was 0.8% of the sum of the normalizations of 
PL1 and PL2. Untying T 3 improved the fit considerably: 
X 2 /dof fell to 9745.1/9159, T 3 was 1.28±g;?|, and A H3 was 
<0.03 xlO 22 cm" 2 (Model 6). All of the improvement in 
the fit when T 3 was thawed came at energies below 2 keV. 
PL3 likely represents a blend of scattered emission (likely 
< 0.8% of the total nuclear hard X-ray continuum) plus 
X-ray emission associated with the kpc-scale jet, with per- 
haps some contribution from XRBs or ULXs in the inner 
2 kpc of the host galaxy. As we have achieved a much bet- 
ter fit with collisional rather than photoionization models, 
we will henceforth assume the soft emission lines are col- 
lisional in nature, although we cannot rule out, e.g., that 
some portion of the O VII and O VIII emission lines may 
be due to photoionization, or that there may be a small 
contribution from a O VII radiative recombination contin- 
uum feature to the observed Fe XVII 3s-2p emission com- 
plex at 0.732t°;°°g keV. For the moment, we adopt Model 
6 as our new baseline broadband model. Best-fit parame- 
ters for Models 5 and 6 are listed in Table 2. Data/model 
residuals for Models 3, 4, and 6 are shown in Figure 5. 

3.1.1. Hard X-ray continuum emission 

Evans et al. (2004) fit the hard X-ray Chandra and 
XMM-Newton data with a dual power-law model, wherein 
Ti and T 2 differed by about 0.3. We untied T 2 from Ti. 
However, compared to Model 6, x 2 om y dropped by 1.4, 
and the best-fit values of Y\ and T 2 were consistent; the 
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uncertainty on F 2 was ±0.24. In subsequent fits, we will 
continue to fix r2=ri. 

There has been no strong evidence from past observa- 
tions for a strong Compton reflection hump in Cen A. We 
added a Compton reflection component to Model 6 using 
PEXRAV, keeping the input normalization equal to that 
of PL 1, the inclination fixed at 30°, the high-energy cutoff 
fixed at 400 keV, assuming solar abundances, and leaving 
the PIN/XIS2 normalization free. The best fit- model pre- 
ferred a value of the reflection fraction R (defined as fl/2n, 
where f2 is the solid angle subtended by the reflector) of 
0, with an upper limit of 0.05, consistent with previous 
suggestions that the Fe Ka line originates in Compton- 
thin material. This limit is identical to that obtained by 
Bcnlloch ct al. (2001) using Rossi X-ray Timing Explorer 
(RXTE) data. Contour plots of R versus Ti and versus the 
PIN/XIS2 relative instrument normalization arc shown in 
Figure 6. 

Using CGRO OSSE, Kinzer et al. (1995) found evidence 
for a high-energy cutoff which appeared at 300 keV to 
~700 keV in the highest- to lowest-flux states, respectively. 
In addition, Steinle et al. (1998) fitted data combined from 
CGRO OSSE, COMPTEL and EGRET; one of the spec- 
tral breaks they found occurred near 150 keV, with the 
spectral index a changing from 1. 74+^06 below the break 
to 2.3±0.1 above it. Keeping the cutoff energies for PL1 
and PL2 in Model 6 equal to each other, we found a lower 
limit of 400 keV to any break in the power law. Future 
improvements to background modeling of the GSO may 
eventually allow us to establish more strict constraints. 

3.2. Hard X-ray emission lines and absorption edges 

In addition to that from Fe, fluorescent Ka lines from 
Si and S have been claimed by Sugizaki ct al. (1997) and 
Evans et al. (2004). We added five more Gaussians, each 
with width tied to that of the Fe Ka line; it was signifi- 
cant at >90% confidence according to the F-test to add 
each line. The centroid energies are consistent with Ka 
emission from Si, S, Ca, Ar, and Ni; detection of the lat- 
ter three lines are reported here for the first time. The 
decrease in \ 2 for each line, including the Fe Ka and K/3 
lines, along with energy centroids, fluxes, and observed 
equivalent widths EW \ )S: are listed in Table 3. In the case 
of Si, we caution that since data near 1.8 keV have been 
ignored due to calibration uncertainties, values of flux and 
equivalent width may be uncertain (e.g., some fraction of 
the large value of A% 2 may be due to calibration uncer- 
tainties). The detection of the Ni line is robust: there 
does exist a Ni Ka line in the NEP background spectrum, 
but it is 200 times fainter than the source spectrum line. 
Data/model residuals to models with all six Ka lines and 
the Fe K/3 line removed are shown in Figures 7-9. 

The Fe Ka line has a best-fit measured width of 14 eV, 
although formally the value is only an upper limit, 24 eV. 
The intrinsic width of the line is then found by subtract- 
ing in quadrature the 55 Fe calibration line width of 9+ 2 eV 
from the measured line width. We infer an intrinsic line 
width of < 23 eV, which corresponds to a FWHM velocity 
of <2500 km s" 1 . 

We next explored the abundances of the hard X-ray ab- 
sorbers (Model 7). The K-shell edges for Fe, Ca, and S 
are detected at >99.999% confidence, determined by set- 



ting each zvphabs abundance value to and refitting. 
The detection of the Ca K edge is likely robust: the re- 
sponse of a line at 6.4 keV includes a Si escape feature 
near 4.5 keV, but it is a factor of more than 50 fainter 
than the continuum and likely does not influence the Ca 
edge. The intrinsically weak Ar K and Ni K edges are 
not detected at high confidence; in the latter case, the de- 
crease in the effective area of the XIS above 8 keV is also 
a factor. We re-fit the data, allowing Zp e , Zc a , and Z$ 
to vary, and keeping the abundances of all three absorbers 
tied to each other for simplicity. We also kept all other 
abundances fixed at solar. Best-fit model parameters for 
Model 7 are listed in Table 2; the edges are also displayed 
in Figures 7-9. The abundance results are listed in Table 
4; note that the Fe abundance is inconsistent with solar at 
3ct confidence. Given the column density of the absorbing 
material and the 2-3 keV rollover, most of the rollover is 
the result of absorption by K-shcll O (with some contri- 
butions to the total opacity from K-shell Ne and L-shcll 
Fe). We are thus actually measuring abundances relative 
primarily to O, not to H. The abundances listed relative 
to H in Table 4 implicitly assume Zq = 1.0. 

To find the equivalent optical depth for each edge, we 
re-fit the data with Zp c , Zc ai and Z§ set to zero in the 
zvphabs components, but with edges at the appropriate 
energies. The results, listed in Table 4, show that the en- 
ergy of each edge is consistent with absorption by neutral 
atoms. 

The Fe K/3/Ka intensity ratio can yield insight into the 
ionization state of the line-emitting gas. However, in the 
present spectrum, the K/3 line is somewhat blended with 
the Fe K edge. However, we note that the Fe K edge energy 
and depth are robust to the properties of the K/3 line, and 
the energy can be used to constrain the ionization state of 
the material, as discussed further in §4. 

We searched for any possible Compton shoulder to the 
Fe Ka line by adding a Gaussian near 6.24 keV, but there 
was no improvement in the fit. We find an upper limit 
of 1.7xl0 -5 ph cm~ 2 s _1 to Compton shoulder emission 
(EW < 10 eV, determined relative to a locally-fit contin- 
uum). 

Grandi et al. (2003) claimed emission due to ionized Fe 
near 6.8 keV using BeppoSAX. We added a narrow Gaus- 
sian at the rest energies of Fe XXV and Fe XXVI, with 
no change in the fit statistic in either case; upper limits to 
the EW of emission lines were 3 and 1 eV, respectively. 

3.3. Abundances of the thermal plasma 

We next explored the abundances Z of the thermal 
plasma relative to solar values. Abundances for the 
VAPEC1 and VAPEC2 components were kept tied to each 
other. We thawed one element at a time from solar, while 
assuming solar abundances for all other elements. Abun- 
dances relative to H can be measured accurately if the 
continuum is well-determined and the thermal emission is 
not too heavily dominated by lines. Given the presence of 
PL3, one must keep in mind that there could be systematic 
effects associated with some derived abundance values. 

We did not find any strong evidence for non-solar abun- 
dances for C, N, or O. The best chance to constrain Zq in 
our model is via the C VI Ly a and Ly (3 lines at 0.368 
and 0.436 keV, respectively. However, there are very few 
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counts at these energies, and the PL3 component domi- 
nates at these energies. Similarly, Zn could, in principle, 
be constrained via the N VII Lya line, but the line is not 
robustly detected. Allowing Zq, to vary resulted in a best- 
fit value of 0.95±0.35 solar, with x 2 decreasing by only 
1.7. We assume Zq, Z^ and Zq, are solar for the rest of 
the paper. 

We did not test the Si abundance: Si XIII emission can- 
not be constrained, as we have ignored data near the Si K 
instrument edge. The best-fit model does predict a mod- 
est Si XIV emission line at 2.05 keV, but at that energy, 
the PL3 component and the absorbed hard X-ray power 
laws dominate. 

In addition to Fe XVII emission lines, Fe is also respon- 
sible for a large bump in the continuum emission from 
about 0.7 to 1.0 keV. Adjusting the Fe abundance causes 
the PL3 component to vary in normalization to compen- 
sate, creating systematic uncertainty to Z-p e derived from 
L-shell emission. We henceforth rely on the abundance 
determined using the K-edge depth (§3.3). 

We did, however, find significant improvements to the fit 
when the Ne or Mg abundances were varied. In the case of 
Ne, the best-fit abundance was 2.71^2 solar; \ 2 decreased 
by 71 (F-test probability 2xl0~ 15 ). For Mg, the best-fit 
abundance was 1.6±0.3 solar, with x 2 decreasing by 19 
(F-test probability 2xl0~ 5 ). 

We re-fit the data with the Ne and Mg abundances in 
VAPEC1 and VAPEC2 thawed (Model 8). We left the 
vapec abundances for C, N, O, Si, and Fe fixed at solar 
(fixing Z-p c at 1.17, the value derived from the Fe K edge, 
had minimal impact). In the best-fit model, Z^ c and Z-^g 
were 2.7±0.3, and 1.9±0.4, respectively. The tempera- 
ture of the VAPEC2 component fell slightly to 0.56±°°* 
keV. With the exception of these parameters, all other 
previously-derived parameters did not change significantly 
compared to Model 7. The best-fit parameters for Model 
8 are listed in Table 2. A contour plot of Nm versus the 
photon index and normalization of the primary power law 
is shown in Figure 10. An unfolded model spectrum is 
plotted in Figure 11. 

3.4. Continuum variability 

Having decomposed the broadband spectrum, it is nat- 
ural to ask whether the observed 2-10 keV variability (see 
Figure 1) is due to variations in the flux of PL1, PL2, or 
both. Figure 12 shows the 2-4 and 5-10 keV light curves, 
summed over all four XlSes and normalized by their re- 
spective means; the 2-4 keV band is dominated by PL1; 
PL3, a component that is not expected to display short- 
term variability, is an average factor of 30 times fainter 
in this band. In the 5-10 keV band, PL1 is a factor of 
roughly 6 times brighter than PL2. Fractional variabil- 
ity amplitudes for these two bands are similar: F var .2-4 = 
3.0 ± 0.2% and F var , 5 _i = 3.5 ± 0.1 %. Moreover, one 
can see that the two mean-normalized light curves display 
virtually identical variability trends. 

The fact that the 2-4 keV band varies implies that PL1 
does vary in flux. If PL2 is constant or only negligibly 

If a lightcurve C consists of a variable component A plus a constant component _B, then F vaT (C) will be equal to -Fvar (A)/(l + where 
fj,A and /ig are the average count rates of lightcurves A and B, respectively. For instance, in the case of Cen A, if PL1 is an average of 6 
times brighter than PL2 in the 5—10 keV band, and PL2 is constant, than -F V ar,5— 10 should be 6/7 X i<var(PLl) (assuming emission from PL3 
is negligible). 



variable, then variability in the 5-10 keV should be di- 
luted relative to that in the 2-4 keV banoP^I. However, 
it is not the case that F vari 5_io is significantly lower than 
F vari2 -4; in fact, F vari5 _io is slightly higher than F val . i2 -4- 
On the other hand, with such small F var values and the 
fact that PL1 is so much brighter than PL2, it is diffi- 
cult to accurately conclude anything about whether PL2 
is constant, as variable, or more variable than PL1. 

While we have confirmed that PL1 varies, the similar- 
ity in F var for both bands and the similarity in variability 
trends means we cannot rule out the idea that both PL1 
and PL2 vary in concert, as expected if a partial-covering 
scenario applies. Further progress can be made through 
continued broadband monitoring spanning a larger flux 
range. 

4. DISCUSSION 

The broad bandpass of the XIS and HXD aboard Suzaku 
has allowed us to obtain a high-quality spectrum spanning 
nearly 3 decades in photon energy. The long exposure 
time, narrow response, and the low background make the 
XIS spectrum one of the highest quality CCD spectra of 
Cen A to date, particularly below 2 keV. 

In §4.1, we compare our results to other recent works to 
test for evidence of long-term variations in several observed 
parameters. The dual power laws that dominate the hard 
X-rays, and the material that absorbs them, are both dis- 
cussed in §4.2. That same absorbing material is likely the 
origin of the fluorescent emission lines; the nature of the 
Fe Ka line is discussed in §4.3. The soft X-ray emission 
from the thermal plasma is discussed in §4.4. Finally, in 
§4.5, element abundances are derived from the hard X-ray 
absorption edges and fluorescent emission lines, the first 
time abundances in the absorbing gas of Cen A have been 
measured. In conjunction with abundances measured us- 
ing the thermal plasma lines, we show that the abundances 
are consistent with enrichment from star formation. 

4.1. Comparison to previous observations 

We can compare our results to previous works to search 
for any variations in flux, photon index, column densities, 
and Fe Ka line flux. 

We find an absorbed 2-10 keV flux of 2.12xl0" 10 erg 
cm~ 2 s _1 , similar to values obtained by Rothschild et al. 
(2006) for RXTE and INTEGRAL observations of Cen A 
between 2003 and 2004. We find a 20-100 keV flux of 
6.4x10" 10 erg cm -2 s -1 , within the range of 20-100 keV 
fluxes measured by Rothschild et al. (2006). We measure 
an unabsorbed 4-7 keV continuum flux of 1.5 xl0~ 10 erg 
cm~ 2 s _1 and a 2-10 keV absorbed luminosity of 2.9x 10 41 
erg s _1 , values almost identical to those found by Evans 
et al. (2004) during the 2001 XMM-Newton observation. 
These facts suggest the source has not varied in hard X-ray 
flux greatly over the last 5 years. 

The photon index, power-law normalization, and ab- 
sorbing column density of the primary (brighter) power 
law, PL1, are all similar to values obtained by Rothschild 
et al. (2006). Evans et al. (2004), modeling XMM-Newton 
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observations in 2001 and 2002 and a Chandra observa- 
tion in 2002, fit partial-covering models (Ti=r2) and dual 
power-law models (T 2 not tied to Ti). In all previous ob- 
servations, the parameters of the primary power law and 
its absorber are very similar to those for the Suzaku obser- 
vation, suggesting a lack of strong variations over the last 
5 years. 

The EW bs of the Fe Ka line in the Suzaku observation 
is consistent with that measured from the XMM-Newton 
observations; we measure the Fe Ka line intensity to be 
2.32 ± 0.1 x 10~ 4 ph cm" 2 s" 1 and the unabsorbed Fe 
Ka flux to be 2.6xl0 -12 erg cm' 2 s _1 , consistent with the 
Chandra and 2001 XMM-Newton observations (Evans et 
al. 2004). Our EW b s values for the Si and S lines are con- 
sistent with those obtained using Chandra HETGS. The 
Suzaku observation thus shows no evidence for variability 
of the Fe, Si, or S Ka lines compared to the 2001-2002 
observations. 

4.2. Continuum emission components 

We now turn our attention to the nature of the two hard 
X-ray power laws, PL1 and PL2, and the less-absorbed, 
soft X-ray power law PL3. The absorbers for PL1 and 
PL2 arc similar in column density to absorbers seen in 
some other radio galaxies (Wozniak et al. 1998). Ccn A's 
dust lane provides an optical extinction of 3-6 mag, which, 
assuming a Galactic gas/dust ratio, corresponds to a col- 
umn TVh ~ 5-10 xlO 21 cm" 2 . The circumnuclear mate- 
rial within a few hundred pc of the black hole, however, 
is believed to have a much higher column, ~T0 23 cm" 2 . 
Both hard X-ray absorbers are thus more likely associated 
with the circumnuclear material than the optical dust lane; 
PL3, meanwhile, is consistent with absorption by the op- 
tical dust lane only. 

Historically, it has not been clear whether the X-ray 
emission is associated with jet or accretion processes. As 
far as jet emission is concerned, several papers in the 1980s 
(e.g., Burns et al. 1983) suggested that the X-ray emis- 
sion was synchrotron emission. However, by constructing 
the broadband SED, Chiaberge et al. (2001) concluded 
that the X-ray continuum was primarily inverse Compton 
emission, with the synchrotron component peaking in the 
far-IR and falling by 10 15 " 16 Hz. We can gauge whether 
one or more of the three power laws observed by Suzaku 
are consistent with X-ray emission from a jet by using the 
radio-X-ray luminosity density correlations of Canosa et 
al. (1999) and Evans et al. (2006), which are based on 
samples of low-redshift radio galaxies. The total 5 GHz 
luminosity density of the nucleus of Cen A is ^6 x 10 20 W 
Hz -1 sr" 1 (e.g., Evans et al. 2006). This value corresponds 
to a 1 keV luminosity density of ~T-2xl0 14 W Hz -1 sr" 1 . 
The 1 keV power-law normalizations for PL1, PL2, and 
PL3 correspond to, respectively, 1 keV flux densities of 
70, 12 and 0.4 yuJy, or 1 keV luminosity densities of ~ 
8xl0 15 , lxlO 15 , and 4xl0 13 W Hz" 1 sr" 1 . It is there- 
fore plausible that PL3 could correspond to jet emission. 
The X-ray luminosity density of PL2 places it an order of 
magnitude above the Canosa et al. (1999) and Evans et 
al. (2006) relations; however, given that there is roughly 
an order of magnitude of scatter in the relations, it is still 
possible that PL2 could be associated with jet emission, 

We can now revisit previous suggestions that there could 



exist two physically distinct sites of hard X-ray contin- 
uum emission. Our best-fit model indicates components 
with unabsorbed 2-10 keV luminosities of 6.7x 10 41 erg s" 1 
(PL1) and 1.3 xlO 41 ergs" 1 (PL2), each independently ab- 
sorbed by different columns. Evans et al. (2004) suggested 
that the primary power law, PL1, arises via Bondi accre- 
tion onto the black hole, consistent with the low accretion 
rate of Cen A (0.2% of Eddington), and consistent with 
PL1 representing X-ray emission far in excess above what 
is expected from jet emission. Evans et al. (2004) also sug- 
gested that their secondary power-law component could 
be associated with X-ray emission from the pc-scale jet. 
Our PL2 component has a 1 keV normalization > 3 times 
higher and an absorbing column density a factor of >10 
higher compared to the secondary power-law component 
detected by Evans et al. (2004). However, a one-to-one 
correspondence between our PL2 and the Chandra/ XMM- 
Newton secondary power law is highly uncertain. We are 
using larger PSFs and extraction radii than Evans et al. 
(2004); the secondary power-law component detected by 
Evans et al. (2004) could be overwhelmed by non-nuclear 
emission (e.g., what we detect as PL3) in the Suzaku spec- 
trum. Meanwhile, Chandra HETGS's low effective area 
> 5 keV makes it relatively insensitive to the presence 
of a power-law component absorbed by a column 7xl0 23 
cm" 2 . It is therefore possible, though not certain, that 
PL2 could correspond to X-ray emission associated with 
the pc-scale jet. In the soft X-rays, PL3 could be the sum 
of diffuse X-ray emission and knots within the kpc scale 
jet, scattered nuclear emission, and emission from other 
point sources resolved by Chandra ACIS. The 0.5-3 keV 
luminosity of PL3 is 1 x 10 40 erg s , similar to the sum 
of the jet and knot component luminosities as reported by 
Kraft et al. (2002) and Kataoka et al. (2006). PL3 may 
additionally reflect uncertainties in modeling the thermal 
gas. 

Alternatively, as suggested by Turner et al. (1997), a 
partial-covering model, with only one site of hard X-ray 
emission, is applicable, as the photon indices of PL1 and 
PL2 are consistent. Our best-fit model is consistent with 
the idea that the central hard X-ray source is powered by 
accretion. 84% of the sky as seen from the central hard X- 
ray source is covered by a column 1.5xl0 23 cm" 2 , and the 
remaining 16% by a column 7xl0 23 cm" 2 . In this model, 
PL3 represents the sum of X-ray emission from both the 
kpc-scale and pc-scale jets, any scattered nuclear emission, 
jet knots, XRBs and other point sources, all modified by 
dust lane absorption, along with any uncertainty in mod- 
cling the thermal gas. 

Finally, we comment on the absence of a strong Comp- 
ton reflection component. The low R value could of course 
indicate the absence of a large amount of Compton-thick 
material that is brightly illuminated by the primary con- 
tinuum and also efficiently radiates the reflection contin- 
uum. An alternate suggestion is that the Compton hump 
is diluted by a power-law component, e.g., from one of the 
jets. However, the large Fe line flux observed is inconsis- 
tent with this notion. In addition, Wozniak et al. (1998) 
demonstrated that spectral fitting of several radio galaxies 
that also displayed weak or non-existent Compton humps 
is inconsistent with this notion. 
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4.3. The origin of the Fe Ka line 

Assuming that the line originates in gas that is in viri- 
alized orbit around the black hole, we can estimate the 
distance r from the black hole to the Fe Ka line-emitting 
gas. The width of the Fe Ka line corresponds to a FWHM 
velocity i>fwhm of <2500 km s _1 . Assuming that the ve- 
locity dispersion is related to wfwhm as <v 2 > — |wp WHM 
(Netzer et al. 1990), assuming a black hole mass Mbh of 
2xl0 8 M Q , and using GA/bh = fv 2 , we find a lower limit 
to r of 6xl0 15 m = 200 light-days. The best-fit edge en- 
ergy, 7.103±0.015 keV, is consistent with absorption by Fe 
I, although ionization stages up to Fe ~V cannot be ruled 
out (Kallman et al. 2004). 

It is plausible that the same material that absorbs the 
hard X-rays is responsible for producing the fluorescent 
lines. The fact that no Compton hump or 6.2 keV Comp- 
ton shoulder are seen, indicating an origin for the Fe Ka 
line in Compton-thin material, supports this notion. On 
the other hand, the material cannot have a column sub- 
stantially less than 10~ 22 cm -2 or else there would be in- 
sufficient optical depth to produce a prominent Fe K line. 
As an estimate of the Fe Ka equivalent width expected in 
this case, we can use the following equation: 

JT P(E)a ph (E)N B dE 

EW ca i c = f c LUf Ka A s — — — (1) 

-r Inline J 

This method assumes an origin in optically-thin gas that 
completely surrounds a single X-ray continuum source and 
is uniform in column density. Emission is assumed to be 
isotropic. Here, f c is the covering fraction, initially as- 
sumed to be 1.0. ui is the fluorescent yield: the value for 
Fe, 0.34, was taken from Kallman et al. (2004). /k q is 
the fraction of photons that go into the Ka line as op- 
posed to the K/3 line; this is 0.89 for Fe I. A is the number 
abundance relative to hydrogen. Initially, we assumed so- 
lar abundances, using Lodders (2003). P{E) is the spec- 
trum of the illuminating continuum at energy E\ En nc is 
the Ka emission line energy. The illuminating continuum 
is assumed to be the sum of both unabsorbed hard X-ray 
power-law components in Cen A, as EW b s values were de- 
termined relative to the total, local continuum. We note 
that the fainter power-law component cannot reproduce 
the entire Fe line as EW b s is too large; the primary power 
law must be responsible for the bulk of the line. Finally, 
Cph(E) is the photoionization cross section assuming ab- 
sorption by K-shell electrons only; all cross sections were 
taken from Veigele (1973^1'). Ah is the column density. 
The two hard X-ray absorbing components in Model 8 
have column densities 1.5 and 7 x 10 23 cm" 2 , which at- 
tenuate power laws with normalizations in a 5.3:1 ratio. 
We therefore use a weighted average Ah column of 2.2 x 
10 23 cm" 2 . (Strictly speaking, the observation in Cen A 
of two different column densities argues against the use of 
a single, uniform column; however, this weighted average 
suffices for our purpose of estimating the total Fe Ka EW 
expected from both absorbers.) 

The value of EW ca \ c is 128 eV, a factor of 1.5 higher than 
EW hs- Assuming that Zp e = 1.17, as the Fe K absorption 
edge depth suggests, yields EW ca i c = 153 eV. Possible ex- 
planations for this discrepancy include the following: (1) 

17 http:/ /www. pa. uky.edu/~verner/photo. html 



The initial assumption of a 100% covering fraction for the 
line-emitting gas may be incorrect. Assuming isotropic 
emission, f c = 0.7 (Z Fc = 1.0) or 0.6 (Z Fo =1.17) may 
fit the data. Such values would be more consistent with 
an origin in a disk as opposed to a shell completely sur- 
rounding the core. (2) The illuminating X-ray continuum 
may be anisotropic. For instance, Rothschild et al. (2006) 
found that Ik was not correlated with Ah and suggested 
that the line-emitting gas may lie along the VLBI jet some 
distance from the black hole. (3) The line-emitting gas 
may be responding to an illuminating flux that is ~0.6 
of the observed flux. This is possible, for instance, if the 
line-emitting gas is >200 light-days from the X-ray con- 
tinuum origin, and does not lie on our line of sight to 
the continuum origin; the observed line flux is a response 
to the continuum flux convolved with a > 200-day time 
delay. Swift/HAT monitoring indicates that the average 
14-195 keV flux of Cen A in early 2005 was very roughly 
15% higher compared to the flux in 2005 August (private 
communcation from the Swift team in 2007(3), possibly 
explaining some of the line flux discrepancy. Continued 
long-term monitoring of the continuum and line fluxes may 
shed light on this issue. Alternatively, the line of sight 
from the continuum origin to the line-emitting gas may 
be blocked, by a clump of gas that absorbs 40% of the >7 
keV flux and that does not lie along our line of sight to the 
continuum. (4) It is plausible that our initial assumption 
of a uniform-column density absorber is incorrect and the 
absorbers are instead clumpy; this is more consistent with 
the observation of two different columns. There could be 
higher columns lying out of our line of sight to the X-ray 
continuum source. For instance, we can use Eq. 1 of Woz- 
niak et al. (1998), which gives the expected Fe Ka intensity 
assuming a cloud with a column > 10 23 cm" 2 lying off the 
line of sight and subtending a fraction Q,/Ait of the sky as 
seen from the source, as a function of the K edge optical 
depth and the spectral index and normalization of the illu- 
minating continuum. We find that a solid angle of Q /4-7T ~ 
0.5 yields a predicted line intensity compatible with what 
Suzaku observes. Such a covering fraction suggests a torus 
or some sort of infinite slab, e.g., a disk- like structure, al- 
though the current data cannot constrain whether such a 
disk would be in the form of a "standard" optically-thick, 
geometrically-thin disk at all radii or a "disk + sphere" 
hybrid (e.g., Esin, McClintock & Narayan 1997), wherein 
the inner parts of a standard thin disk are replaced by a 
radiatively-inefficient flow such as an advection dominated 
accretion flow (Narayan & Yi 1995). 

4.4. Soft X-ray emission from the thermal plasma 

The Suzaku XIS has yielded the best soft X-ray spec- 
trum obtained to date for Cen A, revealing emission lines 
that likely originate in the thermal plasma that extends 
from the nucleus to a radius of ~6 kpc. At least the inner 
parts of the gas are likely heated by the AGN central en- 
gine and/or the nuclear starburst; O'Sullivan, Ponman & 
Collins (2003) demonstrated that the gas in the central 2 
kpc of Cen A has a much higher temperature then the >2 
to 14 kpc gas in the halo. 

Photoionization does not yield as good a fit to the 



See also the SW/t/BAT Transient Monitoring web page, http://swift.gsfc.nasa.gov/docs/swift/results/transicnts/ 
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Suzaku data as collisional ionization models, as the for- 
mer cannot reproduce the strong observed Fe L emission 
complex. Specifically, to obtain a good model fit, we have 
assumed that the emission is entirely collisional in nature, 
and used a two-temperature VAPEC model with k^T 
near 0.2 and 0.6 keV. This is similar to results obtained 
in a wide variety of spiral and starburst galaxies (see ref- 
erences in Strickland et al. 2004). However, we are likely 
sampling thermal emission from plasma spanning a range 
of temperatures. For example, the soft X-rays are likely 
not completely dominated by emission from gas with k^T 
above 0.6 keV, or else we would detect a H-like Mg line 
at least as intense as the He-like Mg line, and the He-like 
Ne line would be extremely faint. Similarly, emission from 
gas with k-gT < 0.20 keV does not dominate, as the Ne X 
would be extremely faint. 

4.5. Constraints on abundances 

We have derived estimates of Zg, Zc a , and Zp c using the 
K shell edge depths, and estimates of Z<j, Z^ e , and ^Mg 
from the thermal plasma (for the remainder of this paper, 
we assume that the abundances in the extended thermal 
plasma and in the circumnuclear absorbing material are 
identical) . 

We can also estimate abundances relative to Fe, e.g., 
Zsi/^Fe, from the Ka fluorescent lines and Eq. 1. Val- 
ues of u> were taken from the X-ray Data Booklel[3- For 
Ni and Ca, /k q = 0.89 and 0.88, respectively, assuming 
neutral atoms (Bambynek et al. 1972). We assume /k« 
= 0.9 for all other elements. The values of EW ca i c , calcu- 
lated assuming solar abundances and a covering fraction 
/ c =1.0, are listed in Table 3. 

Relative abundances are calculated as, e.g., Zsi/Z-p e 
= {EW ohs , si /EW c ^ Si ) I (EW ohs , Fc /EW calc , Fc ), where 
-EWcaic.Fe has been estimated using the abundance derived 
from the Fe K edge depth, Zp c =1.17, and the solar abun- 
dance was used for -EWcaic.Si- We caution that for Ar, Ca, 
and Ni, the lines are weak and the statistics are poor. For 
Si, there may be systematic effects due to current XIS cal- 
ibration as well as from the fact that the continuum near 
1.7 keV is a mixture of all three power-law components. 
Table 5 lists the abundance ratios. There is agreement on 
the values of Zs/Z-p c derived from the edges and from the 
lines. 

We observe [m/H] = log(Zp c ) = +0.1, a slightly higher 
metallicity value than in the metal- rich population of glob- 
ular clusters at 2-20 kpc radii ([m/H] ~ —0.1, Peng et 
al. 2004) or in the outer stellar bulge component (e.g., 
[m/H] ~ -0.2 at a radius of 8 kpc, Harris & Harris 2002). 
There have been previous indications of high metallicities 
in the central regions of Cen A, i.e., indications of over- 
abundances of N and O in young H II regions in the inner 
few kpc (Moellenhoff 1981). Our observed [m/H] value is 
only slightly higher than values for the ISM of other early- 
type galaxies (Humphrey & Buote 2006). It is similar to 
values for the stellar components of large ellipticals in gen- 
eral, [m/H] - 0.0-0.4 (Henry & Worthey 1999). Trager et 
al. (2000) also found the stellar populations of eight field 
ellipticals to have [m/H] ~ 0.0 — 0.4, though it may not be 
straightforward to compare Cen A to field ellipticals due 
to Cen A's recent merger. Our observed value of Zm s /Z-p c 

19 http://xdb.lbl.gov 



is consistent with the value of 2 obtained by Peng et al. 
(2004) for Cen A's globular clusters. 

However, our observed [m/H] value is much higher than 
the average value (of a very wide distribution) found in the 
halo stars of Cen A, [m/H] = -0.4 at radii of 20-30 kpc 
(Harris & Harris 2000, Harris, Harris & Poole 1999). It 
is therefore unlikely that the circumnuclear material and 
the diffuse plasma have directly originated in the relatively 
metal-poor outer halo, unless enrichment via local star for- 
mation has occurred. 

The premerger origin of the circumnuclear material can- 
not be known for certain, although it is plausible that dur- 
ing the merger, gaseous material originally in the outer 
portions of the progenitor galaxies lost sufficient angu- 
lar momentum to be transported to the nuclear regions 
(Toomre & Toomre 1972). Furthermore, connections be- 
tween the nuclear gas and starburst activity at radii of 
hundreds of pc and less are well-known (e.g., Hopkins et 
al. 2006 and references therein). The likelihood that Cen 
A's circumnuclear gas has been enriched due to starburst 
processes is therefore high. 

The degree of metallicity enrichment depends on many 
factors, including the IMF, star formation rate, and the 
total gas consumed, but simulations of disk-disk mergers 
incorporating star formation by Mihos & Hernquist (1994) 
demonstrated that [m/H] values of 0.3-0.4 at the center of 
the resulting merged ellipticals are plausible. The metal- 
licity value derived from the Suzaku data is consistent with 
enrichment due to the ongoing star formation in the inner 
bulge, as a result of the merging process. 

The observed relative element abundances also support 
enrichment. We calculated the expected relative abun- 
dances due to enrichment by Type II and Type la SNe, 
ignoring enrichment via stellar wind processes. We used 
the abundances relative to solar as found in Tsujimoto et 
al. (1995; their Figures 1-2), assuming a Salpeter IMF, 
converting to the abundances of Lodders (2003), and as- 
suming one Type la explosion for every 10 Type II explo- 
sions. For the ratios listed in Table 5, we expect Zo/Zp c 
~2.5, Z-^ e /Z-pe through Zsi/Z Fc to be 1.4-1.7, and Zs/Zf c 
through Zca/Zpe and Z^i/Zp e to all be ~ 0.7-0.9. Most of 
our observed abundance ratios are consistent (or at least 
roughly consistent) with these predictions. The current 
burst of star formation is estimated to have started at 
least 50 Myr ago (e.g., Dufour et al. 1979). Assuming 
star formation at a uniform rate of 30 M yr _1 (Telesco 
1978 measured a star formation rate 10 times that of the 
Milky Way's spiral arms), this is sufficient to enrich the 
gas to the current metallicity levels within that time span 
(M. Loewenstein, priv. comm.), although we cannot know 
for certain how much enrichment of the gas has occurred 
within the last 10 8 yr and how much occurred closer in 
time to the merger. 

Finally, we discuss the argument of Hardcastle, Evans 
and Croston (2007) that "high-excitation" sources, which 
have torus-like circumnuclear absorbing gas and which ac- 
crete at moderate fractions of the Eddington limit (e.g., 
FR lis), are fueled by a supply of cold gas, i.e., gas trans- 
ported to the nucleus as a result of a merger or tidal in- 
teractions. In contrast, "low-excitation" sources, which 
include all FR I's and some of the less powerful FR lis, 
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may be powered by accretion of hot (X-ray emitting) gas. 
As mentioned before, Cen A, although classified as an FR 
I, has heavy X-ray absorption characteristic of most FR 
lis. Our observation of high metallicity, indicative of cold 
gas being transported inwards to the nucleus as a result 
of the merger, is thus consistent with this argument of 
"cold-mode" accretion. 

5. CONCLUSIONS 

The combination of Suzaku's XIS and HXD has allowed 
us to obtain a high-quality spectrum of the nucleus (in- 
ner 2 kpc) of Ccn A spanning 0.3-250 kcV. The long ex- 
posure time, narrow response of the XIS, and low back- 
ground make the XIS spectrum one of the highest qual- 
ity CCD spectra of Cen A to date, particularly below 2 
keV. The HXD PIN above 12 keV is more sensitive than 
RXTE HEXTE or INTEGRAL instruments, allowing us 
to accurately and directly measure the 12-76 keV flux and 
spectral shape. The HXD GSO provides a direct estimate 
of the 45-250 keV flux of Cen A. 

Our best-fit model includes several broadband compo- 
nents. The >3 keV emission can be fit by two absorbed 
power laws. The primary (brighter) power law is absorbed 
by a column of about 1.5x 10 23 cm 2 . The secondary power 
law, with a brightness about 19% of the primary, is ab- 
sorbed by a much higher column, about 7xl0 23 cm~ 2 . In 
our best-fit model, the photon indices were the same, and 
the best-fit value was 1.817^ooio- Untying the photon 
indices did not result in a significant fit improvement. In- 
cluding data down to 0.3 keV, we find that two VAPEC 
components plus a third absorbed power-law component 
can fit the data. The soft X-rays reveal emission lines due 
to He- and H-like O, Fc XVII, He- and H-like Ne, and He- 
like Mg. The strong Fe L emission supports an origin in 
collisionally ionized, rather than photoionized, gas. Two 
VAPEC components of temperatures k^T = 0.2 and 0.6 
keV can model the soft X-ray emission lines well. 

The Suzaku data are consistent with scenarios previ- 
ously forwarded to explain the nature of the various power- 
law components. As suggested by Evans et al. (2004), the 
primary, hard X-ray power law is likely associated with 
Bondi accretion at the black hole. The secondary power 
law could represent X-ray emission associated with the 
pc-scale VLBI jet, although this is not certain. Alterna- 
tively, a partial-covering model, as suggested by Turner et 
al. (1997), may apply: 84% (16%) of the sky as seen from 
the X-ray source is obscured by a column 1.5 (7) xlO 23 
cm 2 . In any event, it more likely that the material obscur- 
ing the hard X-rays is associated with the circumnuclear 
material than with the optical dust lane. The third, soft 
X-ray power-law component is consistent with absorption 
by the optical dust lane only. It likely represents a blend of 
diffuse emission and knot emission from the kpc-scale jet 
(as well as possibly the pc-scale jet), plus emission from a 
population of XRBs and other X-ray point sources, plus a 
likely contribution from scattered nuclear emission, which 
is likely less than 1% of the total hard X-ray power-law 
emission. 

The hard X-ray photon index, the 2-10 keV flux, the 
Fe Ka line intensity and the column density of the ma- 
terial absorbing the primary power law are all consistent 
with recent observations over the last ~5 years, suggesting 



that no large variations in these parameters have occurred 
over the last few years. During the Suzaku observation, 
the hard X-ray flux was not strongly variable on short 
timcscales, showing only a 10% increase in flux through- 
out the observation. 

K-shell absorption edges due to Fe, Ca, and S are signif- 
icantly detected. Several fluorescent emission lines, likely 
originating in the absorbing material, are detected. In ad- 
dition to Fe Ka and Fe Kf3, Ka lines from Si, S, Ar, Ca 
and Ni are detected at at least 90% confidence. The lat- 
ter three are detected for the first time in Cen A. The 
strict upper limit on the amount of Compton reflection, 
R < 0.05, supports the notion that the fluorescent lines 
originate in Compton-thin material. The Fe Ka line width 
is < 2500 km s -1 FWHM. If the line-emitting material is 
in virial orbit, then it is no closer than 200 light-days to 
the black hole. No emission due to He- or H-like Fe is 
detected. The Fe K edge energy is consistent with ab- 
sorption by Fe I; ionization stages above Fe ^V are ruled 
out. The gas does not likely completely surround the X- 
ray continuum source; the equivalent width of the Fe Ka 
line is more consistent with gas which is clumpy and/or 
covers the continuum source only partially, or the bulk of 
the gas may lie off the line of sight to the X-ray source, 
possibly in the form of a disk. 

Metallicities are measured for the circumnuclear absorb- 
ing gas for the first time in Cen A. The high signal-to-noise 
ratio of the XIS spectrum has enabled us to use the K-shcll 
absorption edge depths, relative fluorescent line strengths, 
and thermal plasma emission lines to derive abundances 
(we have assumed that abundances in the thermal plasma 
gas are the same as those for the circumnuclear absorb- 
ing material). We observe a value of [m/H] = log(Zpo) = 
+0.1 for the circumnuclear material/diffuse plasma which 
is much higher than for the halo stars at 20-30 kpc radii. It 
is likely that the merger that created Cen A's famous edge- 
on dust disk also triggered the infall of cold gas into the 
inner regions of Cen A, where it currently accretes onto the 
black hole. Our metallicity observation suggests that the 
accreting material could not have originated in the rela- 
tively metal-poor outer halo unless enrichment due to local 
star formation occurred at some point. The relative ob- 
served abundances are indeed consistent with enrichment 
by a mixture of Type II and Type la supernovae. 

This high signal-to-noise ratio CCD observation has 
yielded breakthroughs in constraining abundances, but 
further progress can be made with a very high spec- 
tral resolution instrument, such as the planned calorime- 
ter aboard Constellation- X. By better constraining the 
widths and intensities of Ka emission lines and soft X- 
ray emission lines, as well as edge energies and depths, 
Constellation- X will thus yield abundance constraints on 
the circumnuclear material in potentially many dozens of 
AGNs once it is launched. 
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Table 1 
Soft X-ray Emission Lines 


Line 
Identification 


Observed Energy 
Centroid (keV) 


Intensity 
(10~ 5 ph cm" 2 s" 1 ) 


EW 
(eV) 


O VII 
O VIII 

Fc XVII 3s-2p ( 3 Pi/ 3 P 2 / 1 Pi) 
Fe XVII 3d-2p { 1 P 1 f s D 1 ) 
Ne IX 
Ne X 
Mg XI 


U.OD4_ 0Q4 

0.654±0.003 

n 7 o 9 +0.004 
u. (csz_ 006 

n Qon+ - 005 
U.82U_ 003 

0.901±0.004 

1.010±0.003 

1 +0.008 
i.ooi_ 004 


10.9 ± 2.1 

9.3+11 
c. o+l.O 
J -°-0.7 
7 4+0-8 
'•*-1.4 
6 S+0.6 

4.1±0.6 
1.5±0.4 


39±0.8 

37+ 4 ° 

23+ 4 

33l 4 6 

36l 3 5 

33±5 

22±6 



Note. — Equivalent widths EW (col. [4]) were determined relative to a locally-fit continuum. 



Table 2 



Best-fit parameters for selected models 



Parameter 


Model 5 


Model 6 


Model 7 


Model 8 


X 2 /dof 


9806.8/9158 


9745.1/9159 


9639.6/9144 


9529.9/9142 


Tl,2 


1 SK9+O.OO8 
i.ooz_ oig 


1 SQl +0.019 
±.oo±_ 010 




1 81 7+0-023 
1 - 81 ' -0.010 


PL1 normalization 1 


n llfi +0.002 
U.11D_ 012 


0.112±0.002 


0.111±0.002 


n , ,,+0.001 
u — 0.002 


A+I,l (10 22 cm 2 ) 


15 2+ ' 1 


15.1±0.2 


14.7±0.2 


14 7+0-3 
i4 -'-0.2 


PL2 normalization 1 


026+ 004 

U.UZD_q 002 


0.022±0.003 


021 + 002 
U ' UZA — 0.003 


0.021±0.003 


N H ,2 (10 22 cm 2 ) 


82±7 


75±8 


72± 1() 


70+^ 


vapecI k B T (keV) 


n ni+0.03 

u - cii -0.02 


0.24±0.02 


0.23±0.02 


n 99 +0.02 
u ' zz -0.01 


VAPECl normalization 


7.2 +o / 1 xl0- 4 


5.0±l\l xlO- 4 


5.0+.^ xlO" 4 


4.3l°'ixl0- 4 


VAPEC2 k B T (keV) 


n 74+0 09 
u -' 4 -0.03 


0.62±0.01 


0.62±0.01 


n 56+ - 01 

u - oo -0.03 


VAPEC2 normalization 


3.7 +( ,J'|xl0- 4 


5.1±0.3xl0~ 4 


5.2+°' 2 xl0- 4 


4.8±0.4xl0~ 4 


r 3 


1 16+ 013 
1 - 1D -0.10 


i 98 +0.08 
J-^O-0.12 


1 41+014 
l - q:L -0.09 


, o, +0.08 
J - ol -0.10 


PL3 normalization 1 


6-9±8"| x !0 -4 


7.4+. / 4 x 10- 4 


7.8±g-| x 10- 4 


7-0±g"| x 10- 4 


JVh.s (10 22 cm 2 ) 


O(fixcd) 


<0.03 


<0.05 


<0.03 



Note. — Best-fit parameters for selected broadband models. All models include two absorbed X-ray power laws (PL1, PL2), 
two VAPEC components, a third soft X-ray power law (PL3), and Fe Ka and K/3 emission lines. Listed are: Model 5, where 
PL3 is not absorbed; Model 6, where PL3 is absorbed; Model 7, where ZVPHABS abundances are thawed for S, Ca, and Fe and 
fluorescent lines for Si, S, Ar, Ca, and Ni are added; and Model 8, our best-fit model, in which the VAPEC abundances Zf^ e = 
2.7±0.3 and Z Ms = 1.9±0.4. 
1 : units are ph keV -1 cm~ 2 s" 1 at 1 keV. 



Table 3 

Fluorescent Emission Line Parameters 





Energy 


Intensity (10 


EW obs 




.F-test 


EW C!llc 


Line 


Centroid (keV) 


ph cm~ 2 s _1 ) 


(eV) 


A X 2 


Prob. 


(eV) 


Si Ka 


1.71±0.01 


0.9±0.2 


24±6 


39.8 


7.6X10" 10 


26.3 


S Ka 


2.307±0.016 


0.8±0.4 


6±2 


5.0 


9.6xl0~ 2 


11.4 


Ar Ka 


2.994±0.023 


1.2±0.5 


13±4 


20.1 


8.0xl0~ 5 


4.7 


Ca Ka 


3.690±0.023 


1.8±0.7 


8±3 


6.2 


5.4xl0~ 2 


3.7 


Fe Ka 


6.394±0.003 


23.2±1.0 


83±3 


3077 


<lxl0- 40 


128 


Fe K/3 


7.056 1 


3.0 2 


ll 2 


28.3 


2.7xl0~ 7 


17 2 


Ni Ka 


7.47±0.05 


1.4±0.8 


15±5 


9.3 


l.lxlO" 2 


7.1 



Note. — Results are for Model 8. Observed equivalent widths (col. [4]) are determined relative to a locally- fit continuum. 
All widths were tied to that for the Fe Ka line, 14^J2 eV. EW CBi \ c (col. [7]) denotes predicted equivalent width values using 
equation 1 (uniform covering assuming line of sight column density), solar abundances, and a covering fraction of 1.0. 

1 denotes a fixed parameter. 

2 denotes value tied to 0.13 that for the Fe Ka line. 
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Table 4 

Hard X-ray Absorption Edges 



K-shell 
Edge 


ZVPHABS 
Abundance 


Edge Energy 
(kcV) 


Edge Optical 
Depth, t 


S 

Ca 
Fe 


1 14+0-11 

i c+0.8 
• L -°-0.6 
-1 17 +0.12 
± - 1 '-0.09 


2.468t°;°it 
4.040±0.045 
7.103±0.015 


24+0- 03 
u -^_0 .06 

0.03 ± 0.01 
0.21 ± 0.01 



Note. — The ZVPHABS abundances in col. (2) are from Model 8. The edge energy and optical depths (cols. [3]-[4]) are from 
a model with ZVPHABS abundances for Fe, Ca and S set to zero and edges fit instead. 



Table 5 
Abundance Ratios 



Method 


Ratio 


Value 


1 


Zo/Zp 


0.8 ± 0.4 


1 


Z^ e /Zp c 


2.3 ± 0.3 


1 




1.6 ± 0.3 


3 


Zsi/Z Fc 


1.5 ± 0.4 


2 




1.0 ± 0.1 


3 


Zs/Zf c 


0.8 ± 0.3 


3 


Zai/Zfc 


4 ± 1 


2 


^Ca/^Fe 


1.3 ± 0.6 


3 


^Ca/^Fo 


3 ± 1 


3 


•^Ni/^Fe 


3 ± 1 



Note. — Methods 1, 2, and 3 (col. [1]) denote abundances derived from thermal plasma emission lines, K shell edge depths, 
and Ka fluorescent emission line intensities, respectively. 
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Fig. 1. — Orbitally binned XIS and HXD light curves. The top panel shows the 2-10 keV count rate light curve, summed over all four XIS 
cameras. The 12-76 keV HXD PIN light curve is shown in the middle panel. The 45-250 keV HXD GSO light curve is shown in the bottom 
panel. 
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Fig. 2. — HXD PIN (left) and GSO (right) spectra. The upper panels show the net source spectrum (gray points), the background (lower 
black points), and the total (source + background) spectrum (upper black points). The PIN spectra have been binned such that the net 
spectrum has a minimum signal-to-noise ratio of IOit per bin. The GSO spectra have been binned as described in §2. The lower panels shows 
the ratio of the net source spectrum to the total spectrum. 
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Fig. 3. — Residuals to Model 1, the single absorbed power law (top, and Model 2, which included dual absorbed power-law components and 
Fe Ka and K/3 line emission (bottom). Black points denote the three FI XIS spectra, which have been co-added here for clarity. Gray, red 
and blue points denote XIS1, PIN, and GSO, respectively. Rest-frame energies are shown. The XIS data have been plotted with a binning 
factor of 10. 



Suzaku Ccn A observation 



17 



1 1 : — i r~ r - r 1 1 1 r— r 




Energy (keV) 

Fig. 4. — Ratio of the soft band data to a simple power law, showing the prominent emission lines. Black points denote the XIS BI; gray 
points denote data from XISO, 2 and 3, which have been co-added for clarity. Rest-frame energies are shown. The data have been plotted 
with a binning factor of 3. 
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Fig. 5. — Residuals to various broadband models are shown. From top to bottom are: Model 3, with soft X-rays fit with a simple power law; 
Model 4, with soft X-rays fit with two VAPEC components only; Model 6, with soft X-rays fit with two VAPEC components plus an absorbed 
power law; Model 7, with K-shell absorption edges and Ka emission lines for S, Si, Ar, Ca and Ni fit; and Model 8, our best-fit model, with 
VAPEC abundances for Ne and Mg thawed. Black points denote the three FI XIS spectra, which have been co-added here for clarity. Gray, 
red and blue points denote XIS1, PIN, and GSO, respectively. Rest-frame energies are shown. The XIS data have been plotted with a binning 
factor of 10. 
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Fig. 7. — Residuals to a simple power-law model in the Fe K bandpass, showing the prominent Fe Ka emission line as well as the Fc K/3 
line at 7.08±0.03 keV, the Ni Ka line at 7.47±0.05 keV, and the Fc K edge. Black points denote the three FI XIS spectra, which have been 
co-added for clarity; gray points denote XIS1. Rest-frame energies are shown. All data have been plotted with a binning factor of 10. 
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Fig. 8. — Residuals to a simple power-law model in the S K and Si K bandpass, showing the Si Ka emission line at 1.71±0.01 keV, the 
S Ka line at 2.307±0.016 keV, and the S K edge. Black points denote the three FI XIS spectra, which have been co-added for clarity; gray 
points denote XIS1. Rest-frame energies are shown. All data have been plotted with a binning factor of 10. 
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Fig. 9. — Residuals to a simple power-law model in the Ar K and Ca K bandpass, showing the Ar Ka emission line at 2.994±0.023 keV, 
the Ca Ka line at 3.690±0.023 keV, and the Ca K edge. Black points denote the three FI XIS spectra, which have been co-added for clarity; 
gray points denote XIS1. Rest-frame energies are shown. All data have been plotted with a binning factor of 10. 
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Fig. 10. — Contour plots of the line of sight column density absorbing the primary power law versus the photon index (left) and 1 keV 
normalization (right) of the primary power-law component in Model 8. Dotted, dashed, and solid lines denote 68, 95.4, and 99.73% confidence 
levels, respectively. The plus signs mark the best-fit values. 
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Fig. 11. — Unfolded broadband spectrum for Model 8, illustrating the three absorbed power-law components and the two VAPEC components. 
All data have been plotted with a binning factor of 10. 
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Fig. 12. — Mean-normalized, orbitally-binned 2-4 keV (black filled circles) and 5—10 keV (gray open circles) light curves, summed over all 
four XlSes. Both light curves display very similar variability trends. 



